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The electronic structure of Cu hexafluoroacetylacetonate, crystallized with a stable nitronyl nitrox- 
ide radical [Ovcharenko et al, Russ. Chem. Bull. 53, 2406 (2004)], is calculated from first principles 
within the density functional theory using the SIESTA method, in two magnetic configurations re- 
flecting parallel or antiparallel setting of S = 1/2 spins of Cu(II) ions to those of organic radicals. 
For a given (high-temperature) crystal structure, the interaction is found to be predominanty anti- 
ferromagnetic, and its magnitude estimated to be 67 cm -1 . This preference is discussed in terms of 
calculated electronic properties (densities of states, molecular orbitals). 



1 Introduction 

Over last decades, a large variety of metalloorganic materials has been synthe- 
sized, many of them having non-trivial magnetic properties. One often denotes 
them generally as "molecular magnets", notwithstanding their affiliation to 
one or another of two large groups (see, e.g., [1,2]). The first is referred to 
as "single molecule magnets" , which are molecular crystals composed of well- 
defined blocks, stable in a solution and remaining well spatially separated (by 
solvent molecules) upon crystallization. Such systems gained much attention 
due to spectacular manifestations of essentially quantum mechanical proper- 
ties (quantum tunneling of magnetization, spin crossiver) in the behaviour of 
ensembles of identical molecules. Another group of systems, although also be- 
ing often refered to as molecular magnets, are in fact polymers, because their 
molecular building blocks gets connected into effectively infinite linear chains, 
or more sophisticated structures. The propeties of a single molecule cannot be 
easily specified in such systems, but a number of other characteristics become 
relevant, e.g. conductivity and long-range magnetic order. Among such poly- 
mers one finds metallorganic ferromagnets with respectable values (>300 K) 
of Curie temperature [3]. 
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Further on, the magnetism is not necessarily due to transition metal (TM) 
ions, but it can be brought about by a free radical as well. The study of poly- 
mers based on stable free radicals has a long history. Mixing up TM ions and 
free radivals in a single polymer system opens possibilities of sophisticated tun- 
ing of magnetic properties. The chemistry offers certain flexibility in bringing 
together the localized spins of TM ions and distributed spins of free radicals in 
a controllable way, in a pattern reproducible throughout the crystal. A recent 
review to this subject has been offered by Blundell and Pratt [4]. 

This poses a question to which extent the electronic and magnetic proper- 
ties of such heterospin systems can be grasped by a first-principles calculation. 
While in general ab initio calculation schemes based on the density functional 
theory (DFT) work well (up to certain limitations) in describing the ground- 
state magnetic structure and certain related properties (magetic anisotropy 
etc., sec Rcf. [5,6] for reviews), heterospin polymer systems have certain obsta- 
cles for a calculation. First, many DFT methods found useful in calculations of 
single molecular magnets, can work on isolated molecular fragments only, and 
do not support periodic boundary conditions. The same applies to more accu- 
rate quantum chemistry calculation schemes. Second, large number of atoms 
and lack of symmetry makes an ab initio calculation with any scheme quite 
demanding. However, a successful calculation may give insight into chemical 
bonding and its relation with magnetic properties in these heterospin com- 
pounds. Moreover one can apply certain constraints in the calculation and 
force the system into one or another magnetic configuration, which would 
help to estimate the typical magnetic interaction energies. 



2 The system of the present study 

We calculate and discuss the electronic structure of one heterospin sub- 
stance, based on copper (II) hexafluoroacetylacetonate Cu(hfac)2, hfac = 
CF3-C(0)-CH-C(0)-CF3, in combination with a stable nitronyl nitroxide 
radical 4,4,5,5,-tetramethyl-2(l-ethyl-l#-pyrazol~4-yl)-imidazoline-3-oxide-l- 
oxyl, described e.g. in Ref. [7]. Some properties of the resulting crystal were 
addressed in Ref. [8]. During last years a large variety of related heterospin 
compounds has been synthesized [7,9,10], in which Cu(hfac)2 blocks are con- 
nected by spin-labeled pirazoles (L Me , L Et , L 1 ~ Pr , L Bu for the above formula 
with different 1-alkyl substitutions). 

The compound in question with the nominal formula C22H21CUF12N4O6, 
refered to as Cu(hfac)2L Et , crystallizes in base-centered monoclinic lattice with 
the space group C2/c (Nr. 15) and lattice parameters a=30.516 A, 6=9.540 A, 
c=25.194 A, /3=123.96°. The primitive cell contains four above given formula 
units (264 atoms). The monoclinic unit cell (hosting twice as many atoms) 
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Figure 1. Monoclinic unit cell of Cu(hfac)2-L. One chain of Cu atoms with their embedding (hfac)2 
fragments and connecting radical molecules runs along the diagonal of the mid-height cutting plane; 
of an (equivalent) chain in the basal plane, only Cu atoms are shown. H and F atoms are omitted. 



is shown in Fig. ^ the arrangement and labelling of atoms in the fragments 
relevant for our discussion, - 'head-to-head' coupling of ligands to magnetic 
Cu atoms, embedded in the (hfac)2 blocks, - in Fig. |2 

3 Calculation method and setup 

The calculation have been done within the DFT, by the Siesta 
method and computer code [11-13], using atom-centered and strictly 
confined numerical functions [14, 15] as a basis set for solving the 
Kohn-Sham equations. Norm-conserving pseudopotentials have been con- 
structed according the Troullier-Martins scheme [16], for the follow- 
ing electron configurations (pseudoization radii in Bohr are given in 
brackets): C2s 2 (1.25)2p 2 (1.25), N2s 2 (1.25)2p 3 (1.25), 02s 2 (1.15)2p 4 (1.15), 
F2s 2 (1.20)2p 5 (1.20), Cu4s L50 (2.08)V(2.08)3d 9 - 50 (1.79). The basis set was 
double-^ with polarization orbitals, according to the notation common in tight- 
binding schemes (see, e.g., Refs. [14,15]) for Cu, N and O, and double-C for C, 
F and H. The energy shift parameter responsible for the localization of basis 
functions in Siesta with its 'default' value of 20 mRy resulted in the following 
maximal extension of basis functions: 2.492 A (H), 2.577 A (C), 2.265 A (N), 
2.083 A (O), 1.899 A (F), 3.308 A (Cu). The sampling over the Brillouin zone 
has been done with 3x5x2 divisions along the reciprocal cell vectors. The 
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Figure 2. An enlarged view of Cu(hfac)2-L-Cu(hfac)2-L sequence, making a period of the chain 

in Fig. with labelling of atoms. 

real-space grid used for the fast Fourier transform in the process of solving the 
Poisson equation was generated with the energy cutoff parameter of beyond 
300 Ry, that resulted in 144x96x240 divisions along the lattice vectors. The 
exchange-correlation was treated in the generalized gradient approximation 
after Perdew-Burke-Ernzerhof [17]. 



4 Results and discussion 

We used the coordinates of atoms as determined experimentally [7, 18] at room 
temperature, without attempting to adjust the structure in the calculation. In 
fact the crystal structure, notably Cul-02 and Cu2-N4 bonds, slightly vary 
with the temperature [7]. Apparently this gives rise to an abrupt change of 
the character of magnetic interaction between Cul and radical spin, from pre- 
dominantly antiferromagnetic (above ~220 K) to predominantly ferromagnetic 
(below this temperature), as the Cul-02 distance grows from 2.237 A (293 
K) to 2.260 (188 K) to 2.281 (115 K). Ovcharenko et al. [7] emphasize that an- 
tiferromagnetic coupling in the axially coordinated >N— 'O— Cu(II)— O — N< 
groups is rather unusual and comes about due to anomalously short Cu-0 
distance at higher temperatures. It was a major objective of our present study 
to verify a tendency for antiferromagnetic cooupling from a first principles 
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Figure 3. Orientation of magnetic moments in two tested AFM configurations. Only Cul atoms 
flanked by O— N- ■ -N— O groups of the radical and Cu2 with their two N4 neighbours are shown. 

calculation. 

The calculation can be initialized with different spin configurations on se- 
lected atoms; this leads to a set of solutions, representing different metastable 
configurations. Two observations helped to reduce the amount of work: i) the 
solutions with opposite setting of magnetic moments of two radicals neighbour- 
ing a certain Cul atom were impossible to bring to convergence, thus indicating 
that a symmetry or antisymmetry of hybridized wavefunctions with respect 
to a central Cul are strongly favoured; ii) the energy difference between solu- 
tions which correspond to various magnetic order between two chains (shown 
in Fig.^), i.e. due to inversion of all spins in one chain, are tiny - at least they 
are not reliably resolvable within the present calculation setup. This under- 
lines an expected fact that magnetic interaction between two (spatially well 
separated) chains are negligibly small. Therefore we imposed antiparallel set- 
ting of spin moments within one chain with respect to the other, and fixed 
the total spin moment of the unit cell to zero, which helped to improve the 
stability and convergency of calculation. This leaves us with several options 
of mutually setting spins of Cul, two radicals coupled to it, and Cu2. We 
keep two most stable solutions, shown in Fig. El and discuss their electronic 
structure in more detail. 

The AFM2 case involves the inversion of magnetic moments of all Cu atoms, 
as compared to AFM1. This leaves the local electronic structure only slightly 
affected; the Mulliken charges and local magnetic moments (Table ^) are the 
same, within 0.01, for AFM1 and AFM2. One finds 1 /iB distributed over the 
radical (referred to as R in the following); similarly, the spin 5=1/2 can be 
associated with each Cu site, even as in reality the magnetic density spills over 
its neighbouring atoms of the (hfac)2 fragment. We found a band dispersion 
over sampled points in the Brilluoin zone to be within 0.1 meV, preserving a 
band gap of ~16 mRy in AFM1 and ~36 meV in the AFM2 configuration. 
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Table 1. Partial charges (Mullikcn populations) Q (e) and local magnetic moments M 
(fJ>&) in two most stable magnetic configurations. Four oxygen atoms in each Cu(hfac)2 
fragment, labeled 03 and 04 in Fig. [2] are in reality subdivided into two inequivalcnt 



sites, according 


; to inversion 


symmetry relative to Cu 


centers. 








Cul 


04 


Cu2 03 


Ol 


Nl 


N2 


02 






AFM1 










Q: 10.64 


6.15 


10.67 6.14-6.16 


6.10 


4.93 


4.90 


6.15 


M: 0.51 


0.07 


0.57 0.08-0.09 


0.29 


0.21 


0.20 


0.19 






AFM2 










Q: 10.65 


6.15-6.16 


10.67 6.15-6.16 


6.10 


4.93 


4.90 


6.15 


M: 0.45 


0.08-0.04 


0.57 0.07-0.08 


0.30 


0.21 


0.21 


0.18 



Although not a strict criterion, a larger band gap is often an indication of a 
more stable configuration among competitive ones. Indeed, the total energy is 
lower in the AFM2 case by 35 meV, per unit cell of 256 atoms, which contains 
two Cu2-L-Cul-L units. It could be expected that the experimental band 
gap be larger, because the molecular orbitals just below and above the gap 
are those with large Cu contribution, and it is known that an inclusion of 
correlation effects within the (quite localized, in this case) Cu3d shell beyond 
the conventional GGA would tend to move these states apart. 

Assuming the Cul-R magnetic interaction as dominant and those involv- 
ing Cu2 as comparably small, one can map the above energy difference onto 
the Heisenberg model. There are several ways to extract magnetic interac- 
tion parameters from a first-principles calculation, e.g., in a form consistent 
with spin-fluctuation theories, in terms of Green's function elements [19-21], 
or via non-local susceptibility [22]. The difficulty in the present case would 
be a delocalized character of a spin over four atoms (Ol-Nl- • -N2-02) of the 
radical, and a certain ambiguity in defining the corresponding matrix elements 
of the Green's function and/or susceptibility. A less conceptually problematic, 
although quite 'old-fashioned', is comparison of total energies in several rep- 
resentative magnetic configurations. A prerequisite for this approach to be 
meaningful is that the involved spins retain their 'identity' (i.e., magnitude, 
localization etc.) in the configurations to be compared. This is indeed the case 
in the compound under study, as follows from the inspection of local densities 
of states (DOS, shown below) and spatial spin densities. 

Given ambiguities in the definitions of sign and prefacor of the Heisenberg 
model, we introduce the magnetic interaction parameter consistently with 
Ref. [10], notably Eq. (2) therein for the special case of interaction within 
the >N-'0-Cu(II)-0'-N< 'exchange cluster': 

TL = —2JScuS' + (Zeeman term) + (interaction with more distant ions) , 

where S' is the total spin of two radicals flanking Cul. Mapping this onto 
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total energy values of our calculation yields: 

-Eafmi — -Eafm2 = — 16JSc u Sr , 

whence (with Sr=1/2 for the spin of one radical) —J = 8.3 meV = 96 K = 
67 cm -1 . This rough estimate of (antiferromagnetic) coupling between a Cu 
ion and the radical can be looked at as an upper-bonded value: an inclusion of 
intraatomic correlation effect beyond the conventional GGA may scale down 
the energy difference between magnetic configurations by a factor as large 
as 3-4, depending on a system, according to a previous experience for other 
molecular magnets [23,24]. Moreover, an explicit inclusion of interactions with 
Cu2 into the model would slightly renormalize the Cul-R interaction further 
down. 

We could not find an experimental value to compare with, e.g. that ex- 
tracted from a Heisenberg-model fit to magnetic measurements data; the value 
J=26 cm" 1 given in Table 3 of Ref. [7] stems from the fit over low-temperature 
range, where exchange interaction within the R-Cul-R group is of predomi- 
nantly ferromagnetic character. (We remind that our calculation has been done 
in the room-temperature structure, whereas at 220 K the exchange interaction 
does apparently change sign). Antiferromagnetic type of coupling over large 
temperature range was observed in chemically related, but structurally differ- 
ent, system Cu(hfac)2L Pr , yielding J=— 100 cm" 1 (see Ref. [7] for details). It 
is noteworthy that all estimations of "additional" exchange parameters (those 
incorporating, on the average, interactions to more distant Cu2 atoms and 
between the chains) are by almost two orders of magnitude smaller. This is 
consistent with relative unimportance of these more distant interactions also 
according to the calculation. 

An inspection of partial DOS helps to understand the origin of energy pref- 
erence of the AFM2 configuration. Fig. 0] compares the partial DOS at Cul 
and its neighbours which exhibit a noticeable magnetization, and that of more 
distant Cu2. For better comprehensibility, a very structured DOS from nondis- 
perive energy bands was broadened with half-width parameter of 0.1 eV. Ob- 
viously this completely smears the above-mentioned tiny band gap. Although 
formally dealing with a non-metal, we refer to the energy which separates 
occupied and unoccupied states as the Fermi energy (zero in Fig. 

Somehow simplifying, but in accord with chemical intuition, the Cu(II) ion 
has its 3d majority-states fully occupied, and in the minority-spin one state 
vacant. This vacant 3d state is responsible for a narrow peak in the DOS 'at' (in 
fact, just above) the Fermi level. This simple picture seems not quite consistent 
with nominal magnetic moments at Cul and Cu2 sites as shown in Tabled 
but that's merely because the Cu magnetization spills over its neighbours. It 
is well seen in spin density maps (not shown here) and in the maps of Cu- 
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Figure 4. Spin-resolved partial DOS in two AFM configurations. Zero energy separates occupied 
and unoccupied states. The features marked by thick dots and labeled a and b are explained in the 

text. 



centered molecular orbitals shown below; moreover it is demonstrated by a 
presence of the Cu-related peak just above the Fermi level in the 04 DOS. 
Apart from this, we have a magnetic moment of the radical, roughly coming 
from summing up the values in Tablen° ver 01, Nl, N2 and 02 atoms where 
this moment is distributed (with a pronouncedly larger portion residing on the 
more Cu-distant 01). 

As a most coarse observation, the Cul-DOS is 'inverted' when changing 
from AFM1 to AFM2, but this is accompanied by a noticeable reconstruction 
of electronic states, responsible to the Cu3d-02p hybridization. Among Cu- 
related and R-related orbitals, which share the energy interval of —4 to +1 
eV around the Fermi level, we further concentrate on some, with particular 
character of bonding. The highest-energy peak among the fully occupied ones 
in majority-spin panel on the Cul site (spin-up for AFM1, spin-down for 
AFM2) comes from a state which is nonbonding with respect to the radical and 
antibonding with respect to oxygens of the Cu(hfac)2 group (note its presence 
in the 04 spin-up DOS). The inversion of Cul spin changes merely the position 
of this peak: from -0.83 eV in AFM1 to -0.68 eV in AFM2. Further on, there 
is a state nonbonding with respect to Cul and pinned at the radical. It emerges 
in the DOS of 01, Nl, N2, 02 as the lower of two peaks crossed by the Fermi 
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level. Its wavefunction is shown in Figs. [5] and as that of the last-but-highest 
occupied molecular orbital (-HOMO). It contains no contribution from Cul 
whatsoever (having in fact a node at the Cul position) and is almost identical 
for AFM1 and AFM2 cases. Then there are two pronounced states which 
involve both Cul and the radical: (a) a feature in the minority-spin DOS at 
the Cul site, i.e., in the spin-down channel of AFM1 (at —0.71 eV) and in 
the spin-up channel of AFM2 (at —1.13 eV); (b) the 'just occupied' part of 
the highest of those two peaks which make a distinguished two-peak structure 
best seen in the local DOS of the radical atoms, crossed by the Fermi level. 
In fact this is the HOMO. This state is antibonding with respect to Cul-02 
interaction. These (a) and (b) features are marked by thick dots in the DOS 
figures. 

The lowest unoccupied molecular orbital (LUMO) centered at the selected 
Cul site corresponds to spin-down channel (the spin-up LUMO resides at a 
Cul atom in the other polymer chain) and contains no contribution from R. 
This is understandable, because within ±0.5 eV around the Fermi level there 
is a gap in the spectrum of spin-down states of the radical. The structure of 
the Cu3o? wavefunction within the LUMO is of the xy type, in contrast to 
spin-up z 2 in the HOMO. The next to LUMO (and higher by 16 meV) is the 
(also spin-down) orbital centered at Cu2. 

When we now turn to the molecular orbitals of the AFM2 configuration 
(Fig. H)J) we notice that the last-but-HOMO orbital, nonbonding to Cul and 
antisymmetric relative to its position, has the same shape as for AFM1. How- 
ever, the HOMO differs from that of the AFM1 case by its symmetry at the 
Cul site. It bounds two 04 atoms of the Cu(hfac))2 group much stronger than 
two others; this contributes to much higher disproportionality of magnetic mo- 
ments over 04 atoms in AFM2 case than in AFM1 (see Table [TJ. The LUMO 
of AFM2 does not reflect the Cul states nor those of the radical, but instead 
the Cu2 (see below). The next-to-LUMO contains again the Cul contribution 
in somehow similar shape as in AFMl-LUMOj, spread over the Cu(hfac)2 
plane. However - and this is crucial - the discussed wavefunction spreads over 
the 02-N2- • Nl-01 group, in antibonding setting relative to Cul. 

As a result we see markedly stronger interaction of Cu3d states with the 
magnetic states of the radical in the AFM2 configuration than in AFM1. 
This is so because the responsible states which tend to interact, based on 
their spatial shape, belong now to the same spin channel, namely to spin-up 
in the notation of Fig. 0J As a single manifestation of this effect, the Cu-R 
hybridized state (a), on the one side, and (b), but also the '(LUMO+l)f of 
Fig. El on the other side, move apart. A similar argument applies with respect 
to other hybridized Cul-R states, resulting in a slight overall (albeit non- 
uniform) replacement of occupied bands downwards. As a result, the AFM2 
configuration wins over AFM1 due to the gain in band energy; moreover, this 
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AFM1 




(HOMO-l)t HOMOT LUMO| 

Figure 5. Highest occupied and lowest unoccupied molecular orbitals corresponding to fc=0 in the 

AFM1 configuration. Only the Cu(hfac)2 block and its two neighbouring ligands are shown. 
Positive and negative values are indicated by two grayscale levels. As in an AFM configuration all 
orbitals are doubly degenerate, and the figures show the vicinity of a selected Cul atom, the spin 
label is given for those respective orbitals which give a major contribution at the site in question, 
consistently with Fig. [3] The opposite-spin orbitals at the same energy are prominent at different 
Cul atom in the other polymer chain (see Fig- El • T ne plots are created with the XCrySDen 

program [25]. 

accounts for a slightly larger band gap in AFM2. 

The Cu2-DOS (the bottom panels in Fig. |1] seems largely unaffected by a 
reversal of its magnetic moment: its 3<Z-j-n) DOS in AFM2 is almost identical 
to 3<ij(j) one in AFM1. There are tiny differences responsible for the coupling 
of Cu2 to the 'tail' of the ligand via N4 (see Fig. [2]), but they are left beyond 
the scope of the present study. It suffices to mention that the interaction with 
differently magnetized Cu2 changes the sequence of orbitals in the vicinity of 
the Fermi level, which are anyway nearly degenerate. In the AFM1 configu- 
ration, the LUMO has a pronounced contribution at the Cul site (shown in 
the right panel of Fig. EJ), with the spin-down orientation, and the next state 
(higher by 16 meV) is that centered at Cu2. In the AFM2 configuration, the 
state with Cu2 predominance becomes the LUMO and have no contribution 
at Cul, whereas the Cul-centered orbital is now the next after LUMO (by 4 
meV higher, shown in the right panel of Fig. EJ), now evidently in the spin-up 
channel. 

Summarizing, we investigated electronic structure, local densities of states 
and the shape of molecular orbitals responsible for chemical bonding and mag- 
netic interaction in a heterospin polymer, connecting Cu(II) and nitronyl ni- 
troxide radical in a 'head to head' configuration. We reproduced in the calcu- 
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AFM2 




(HOMO-l)t HOMOT (LUMO+l)T 

Figure 6. The same as in Fig. for the AFM2 configuration. 



lation and explained at microscopic level the preference for antiferromagnetic 
coupling between Cu and radical spins, and made an estimate of the magnitude 
of this coupling. There are certain issues which seem appealing for subsequent 
study. First, one should elucidate the interaction between Cu2 ion coupled 
to the 'tail' of the ligand and the presently analyzed 'exchange cluster'. Sec- 
ond, the variation of all interaction parameters with temperature, according 
to (available by now) crystallographic data [26]. This might help to under- 
stand an unusual switch in the character of magnetic interaction, apparently 
accompanying a structural transition at 220 K described in Ref. [8,27]. 
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